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Abstract 



We present the solution of the time development of an unstable initial density perturbation in 
the linearized Vlasov equation, completing the previous analysis in the literature. The additional 
contributions found are usually damped and can be neglected at large times in the unstable region. 
£Sj ■ The work clarifies also the problem of the normalization of the solution with respect to the initial 

^ ' perturbation of the density. 
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Q\ • With the development of the stochastic transport equations for nuclear collisions , i.e. the Boltzmann- 

Langevin equation, the description of the multi fragmentation in intermediate energy nuclear collisions 
has been challenged 0, ||, The spinodal decomposition mechanism of fragment formation is based 
on the idea first proposed by Heiselberg et al. Q] that the nuclear matter at low density p ~ Po/3 
is unstable against density perturbation. In the linear regime, any perturbation with unstable spatial 
momentum k will develop in time as e r ( fc )*. The fastest growing mode will dominate the linear response 
and consequently, the typical fragments will have the size 1/fc , determined by the largest growth rate 
r(fc). The importance of the nonlinearities of the Vlasov equations for the fragment formation is still not 
settled ||. However, even the linear response analysis presented so far are often incomplete. 

It may seem that after many works dealing with the development of the instabilities in the Fermi 
systems everything is understood in the linear regime [| |[ [l| . Usually, people study the proper modes 
of the Fourier transformed linearized Vlasov equation : 

i(kv - to)Sf{uj, k,p) - i^^kv^-8p(u, k) = , (1) 

where Sf(cu, k,p) — J dt d 3 x 6f(t, x,p) e Mt ~ jkx is the Fourier transformed in time and space perturbation 
of the equilibrium phase-space distribution no, Uk is the Fourier transformed density dependent mean 
field. In the spinodal region the eigenvalue problem for the linearized Vlasov equation has two imaginary 
frequency solutions u>(k) — ±iT(k) in certain range of k < k max limited by the range of the mean field 
potential. The solution using the eigen-functions of the Vlasov equation is not a well posed problem for 
the solution of the time development of the initial perturbation (or for the solution of the Vlasov equation 
with a noise term Q). The analysis using one sided Fourier transform are more suitable and are a well 
posed boundary value for the Vlasov equation. In la], using the one-sided Fourier transform only the 
large time limit is explicitly discussed. References jj], [|] discuss the singularities corresponding to the non- 
collective modes. In the case of discrete momentum these singularities correspond to the cut singularity 
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in the continuum momentum studied below. In the framework studied in this letter the calculations of 
the contribution from the non-collective modes can be easily done. In particular we can show that in 
non-exceptional cases of the initial perturbation the additional contribution to the density development 
are damped. In the discrete formulation j^] it requires the use of a large basis of as complete as possible, 
since the damping can arise only from the beating in frequency of the sum of many contributions with 
real frequency. 

Below, we shall calculate the evolution of the initial density perturbation Sf(t = 0,k,p) = g(k,p) in 
the linearized Vlasov equation in the unstable region. Taking the one-sided Fourier transform : 

POO 

6f(ui,k,p)= / dtf(t,k,p)e iut , (2) 
Jo 

one obtains : ^ 

z(kv - u))Sf(u, k,p) - i——kv—^-6p(uj, k) = g(k,p) . (3) 
op Oe 

Dividing by kv — to and integrating over p an equation for the density perturbation is obtained. The time 
dependence of the density perturbation Sp can be found using the inverse transform : 



2n e(k,Lu) 



with 



G(k,oj) = / T-—^- , (5) 



(2tt) 3 kv - 



ui 



and 

e(fc,w) = l-— /-—,———— , (6) 



dp J (2n) 3 kv — lo de 

the integration path in the inverse Fourier transform laying above any singularity of the integrand. The 
integral can be calculated closing the integration path in the lower half-plane. The solutions of the 
dispersion relation : 

e(fc,w)|* = (7) 

gives the poles of the integrand. The contribution to the integral from the residues at these poles gives 
the eigenfunction solution of the linearized Vlasov equation. In the spinodal region, for the wave-vector 
k of an unstable mode the dispersion relation (0) will have two solutions with imaginary frequencies. 
The two frequencies give a damped and a growing component in the time dependence of the density 
fluctuation : 

An d h\ - V" ~G{k,±iT{k)) ±r(fc)t 

dp po ie(t,k) - > 7r-7T — r"7Tj j e . (8) 

For large times the part with growing exponent dominates. However, this solution is incomplete. For 
t=0 we do not recover the initial condition : 

f d 3 p 

Sp po i e (t = 0,k)^ J J^^9(k,p) , (9) 

e.g. for 

g(k lP ) = (2ir) 3 g(k)5 3 (p-po) (10) 

at zero temperature one gets for the overlap of the density disturbance with the eigen-functions of the 
Vlasov equation : 

Sp po ie(t = 0, fc) -2F(fc) 2 (r(fc) 2 + {kvff) 

g (k) (r(fc)2 + (fc W/ )2(i + F (fc)))(r(fc)2 + (kp )2) ? > y > 
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where Fo(k) = ^-^f is the Landau parameter Fq of the Landau- Fermi liquid, pf and Vf denote the 
Fermi momentum and Fermi velocity respectively. Also generally, the sum of the overlap with the stable 
and with unstable eigen-functions is different from 1 . The original integral in the right hand of (^) has 
a correct limit at t = : 

00+4,7 diu G(k,u) _ f d 3 P 
OQ+iCT 2tt e(k,w) ~7 (2tt) 3 



Sp(t = 0,k) = -i I —-£±-L=l*g( k ,p), (12) 



as one can convince oneself calculating the residue of the integrand at infinity. The discrepancy between 
the two results arises from the contribution of the cut in the integrand on the right hand side of (^) on 
the real axis in frequency w. 

The contribution from the poles of the inverse susceptibility X/e(k,u>) to the time development of the 
density is given in eq. (|J). However, the two poles are not the only singularities of the inverse suscepti- 
bility. The singularity structure of the susceptibility on the real axis can be read off after performing the 
angular integration in (^]) : 

e(k,u) = l-F {k) f°° x dx (13) 
J xkvf ox 

where $(s) = 1 — ^ n (j~r) an( i where in the rescaled variable x the equilibrium distribution is fo(x) — 
exp(e F (x'^-i)/T)+i ' ^ ne logarithm in <fr(s) in the the susceptibility has a cut on real axis : 

e(fc, w±ie) = Re e(k,ui) ± ilm e(k,u) 

dfo(x) 



F (k) / x dx Re<$>(^—)- 

Jo y xkv f 



dx 



Also the function G(k,u>) can only have singularities on the real axis. Accordingly, the contribution of 
the cut to the density perturbation is : 

*-m) - - r - 1^>-" . da) 

J_ 00 2n \ e(k,uj + le) e{k,w - le) J 

This is different from the cut singularity in relativistic Vlasov equation where the integration region is 
always finite. 

Taking the perturbation well localized in momentum ( |l0| ) and using kvo ^ UJTi(i — P kv '_ a ±iTr5(u>— kvo) 
we obtain : 



6p(t,k) = —pf 



ds Im e(k,skv f ) e -iakv f t 

so — s (Re e(k, skvf)) 2 + (Im e(k, skvf)) 2 



(Re e(k, s kv f )) 2 + (Im e(fc, s kv f )) 2 ' V ' 

where so = j^j . Except for the k dependence of the Landau parameter Fo(k), the time scales appearing 
in the time dependence of the density Sp cu t(t,k) are proportional to 1/k. It is the same scaling as for 
the growth rates T(k) ~ k for constant Fq . In Fig. 1 we show the time dependence of <5pc ^*' fc - > at T = 
and at T — e//2. The evolution of the density dp cu t(t,k) as given by the cut contribution is oscillating 
and undamped for initial perturbation with definite momentum. The phenomenon is similar to the Van 
Kampen modes in the relativistic Vlasov equation 0, |f| . Apparently, these modes have the overlap 
picked only on a finite number of eigenfunction in the expansion in ref (2| . 

Generally, the initial density disturbance will be a function of the momentum g(k,p). Although such 
a function can be written as a superposition of 5 functions in momentum, its behavior in time will be 



3 



different. The phenomenon is analogous the beating in frequency mentioned above. Let us study the 
evolution of a deformation of the Fermi surface : 

f/(t = o,M=*(*)^^. d7) 

PfK oe 

where Kq — fo(x)dx. The time development of the density disturbance is given by the sum of the 
term from the poles of l/e(fc, w) and form the cut on the real axis in frequency. The projection on the 
eigen-functions of the density G(±iT(k),k) will grow or will be damped exponentially according to eq. 
(||). The overlap of the initial density perturbation with the eigen-functions of the Vlasov equation as 
a function of the instability rate is shown in Fig. 2. The overlap is greater than 1 with the limits 1 and 
2 in the case of strong and weak instability respectively. This shows that the contribution from the cut 
cannot be neglected at small times in order to insure the correct limit at t = 0. The contribution from 
the cut can be found from eq. ([lq), with : 



G(u>±ie, k) = ~ rtfs } — / x dx In 



2K kv f 



LU + xkvf 



LO — xkvf 



dfo(x) 



dx 



j*f**. us) 

2K kv f J\ u \/ kVf dx 

The result of the numerical integration is shown in Fig. 3. We see that the contribution form the cut 
to the density distribution is rapidly damped in time. As in the previous cases the time scales involved 
in the time dependence of the density disturbance are proportional to the wave-vector fc, for constant 
F . For both cases of the initial distribution studied we obtain the correct limit of the total density 
distribution at t = 0, 5p po i e (t — 0, k) + 5p cut (t = 0, k) = g(k) . The finite temperature the curve is 
smooth, almost exponential, without the oscillations found at zero temperature. These zero temperature 



oscillations reflect the sharp limits of the integration in the formula (15), since at zero temperature the 
cut on the real axis is limited to the interval \lo\ < kvf. 

We have found that the description of the time development of an initially unstable density perturba- 
tion can be described as a sum of two terms. The term coming from the poles of the inverse susceptibility 
function, with a growing and a damped mode. Unlike in the stable region the discontinuity in the suscep- 
tibility function does not modify the result for this contribution to the density evolution. In particular 
the position of the zeros of the susceptibility function is not changed, in contrast to the appearance of the 
Landau damping for the stable mode of the Vlasov equation jj^j . A second contribution was found coming 
from the cut on the real axis in frequency for the susceptibility function e(k,u>) and for the one-sided 
Fourier transform of the density perturbation G{k,Lo). We have shown that this contribution is bounded 
in time, and generally strongly damped. Only in the case of a singular perturbation in momentum an 
oscillating solution for the contribution from the cut was found. This contribution can be than neglected 
at large times, where the unstable mode dominates. However, as it can be seen in Fig. 4, it is important 
at small times. In particular the contribution from the cut integration to the density disturbance is 
required in order to recover the initial density perturbation from the solution of the Vlasov (U) equation 
at zero time. Indeed, the overlap of the initial perturbation with the growing and damped mode is mostly 
larger than 1 . It is compensated by a negative contribution to the density at zero time coming from the 
cut integration. It would be incorrect to normalize the pole contribution to the density disturbance at 
zero times Sp po i e (t — 0,k) to g(k) . The correct limit at t = of the time development of the density 
perturbation originates from the following sum rule : 

f d 3 P ^ -G(k,±iT(k)) 



if 



de(k,uj)/duj\ ul=±ir{k) 

dui ( G(fc, lu + ie) G(fc, to — ie) 
_ QO 2tt \ e(k, lo + ie) e(k, to — ie) 

d 3 P f(p) 
(2tt) 3 kv - to 



(19) 



G(M)=/7^^. (20) 
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This sum rule corresponds to taking the complete basis for the expansion of the initial perturbation in 
ref§. 

The results presented in this paper, explain how the initial density disturbance gets additional growing 
term corresponding to the rapid decay of the negative contribution from the cut integration. This effect is 
the most pronounced at small values of the growth rate T/k, where the overlap of the density perturbation 
with the eigen-functions is larger than 1, although finite, ft modifies the result of a simple growth of 
unstable eigen-functions of the Vlasov equation by the value of the overlap at the initial time. This 
enhances the initial perturbation for the wave- vectors corresponding to small growth rate T/k. It means 
that that the prediction that the most unstable modes would be preferably excited would be weakened, 
since the overlap with unstable eigenfunction is close to 1/2 for large T/k. The contribution form the cut 
integration can be neglected in this case. On the other hand, it gives a negative and strongly damped 
contribution at small T/k, which results in a faster increase of the density disturbance at small times. 
At large times, this effect is taken into account as a normalization of the overlap of the initial density 
perturbation with the unstable mode different form 1/2 . 

We have presented the formula and figures in a general way, using variables scaled by 1/kvf. It shows 
that, except for the k dependence of the Landau parameter Fq, the time scales of the oscillation or decay 
of the contribution to the density disturbance coming form the cut integration scale with 1/kvf. Similarly 
as the growth rate T(k) ~ kvf. This means that the decay times of Sp cut go to infinity as k — ► 0. 
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Fig.l 

The time development of the initial density perturbation singular in momentum (10) (Fq = —1.5, 
kvo = O.bkvf). The solid line and the dashed line represent respectively the real and imaginary part 
of the density perturbation at zero temperature. The dotted line and the dashed-dotted line represent 
respectively the real and imaginary part of the density perturbation at T — e//2. 




Fig. 2 

The sum of the overlap of the density perturbation ( |l7| ) with the eigen-modes of the Vlasov equation 
as function of the growth (damping) rate T at zero temperature. Very little change occurs at non-zero 
temperature. 
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Fig. 3 

The time development of the initial density perturbation(|l7|) (Fq = —1.5). The solid line and the dashed 
line represent the result at zero temperature and at T = e//4 respectively. 
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Fig. 4 

The time development of the initial density perturbation ( |l7| ) at zero temperature (i<b = —1.5). The 
solid line represents the pole contribution. The dashed line represents total density perturbation, i.e. the 
sum of the pole and the cut contribution. 
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